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Abstract 

In a continuous-wave terahertz system based on photomixing, the measured ampHtude of the 
terahertz signal shows an uncertainty due to drifts of the responsivities of the photomixers and of 
the optical power illuminating the photomixers. We report on a simple method to substantially 
reduce this uncertainty. By normalizing the amplitude to the DC photocurrents in both the 
transmitter and receiver photomixers, we achieve a significant increase of the stability. If, e.g., the 
optical power of one laser is reduced by 10%, the normalized signal is expected to change by only 
0.3%, i.e., less than the typical uncertainty due to short-term fluctuations. This stabilization can 
be particularly valuable for terahertz applications in non-ideal environmental conditions outside of 
a temperature-stabilized laboratory. 

PACS numbers: 300.6495, 120.6200, 040.2235, 230.0250 
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I. INTRODUCTION 



Terahertz spectroscopy requires a reliable determination of the amplitude and phase of 
the terahertz radiation. One established method for the generation and coherent detection 
of continuous-wave terahertz radiation is photomixing [l-7|. Typically, terahertz radiation 
is generated by illuminating a biased photomixer (the transmitter, Tx) with the optical 
beat of two near-infrared lasers, and coherent detection is achieved by measuring the pho- 
tocurrent Jthz in a second photomixer (the receiver, Rx) which is illuminated by both the 
optical beat and the terahertz radiation. The accuracy of the quantity of interest, e.g. the 
transmission Tgampio of a given sample, depends on the stability of the measured terahertz 
signal /thz- This signal in turn depends on the stabilities of the responsivities S^x and S'rx 
of the two photomixers and on the stabilities of the optical powers of the two lasers at the 
two photomixers. All six of these quantities are sensitive to temperature drifts. Here, we 
report on an efficient normalization of the terahertz signal by employing both photomixers 
as powermeters. 

II. CONCEPT 

Our aim is to minimize the uncertainty of /thz- In order to motivate an expression 
describing /thz, we start from the generation of the terahertz wave at the transmitter. For 
a constant bias voltage of f/bias ^ 10 V, the transmitter on the one hand behaves like a 
photoconductive resistor, i.e. the DC photocurrent at the transmitter, Jtx, is given by the 
bias-dependent responsivity S'Tx(^bias) and the total incident optical power Ptx = -Pl,Tx + 

-f2,Tx, 

/tx = 5Tx(f/bias) X (Pi,Tx + P2,Tx) , (l) 

where Pi,tx and P2,Tx denote the optical powers of the two lasers at the transmitter. On the 
other hand, the interference between the two laser beams gives rise to a beat which yields a 



hus the transmitter emits a terahertz wave 
. In transmission geometry, the power at 



current oscillating at the difference frequency, 
with power Pthz.Tx oc S'|^(f/bias) x Pi,tx-P2,Tx 
the receiver is given by Pthz,Rx oc Tgampie x Pthz.Tx- The photocurrent in the receiver, 
Ithz = -^o.THz cos(A(y9), depends on the phase difference Aip between the terahertz wave 
and the optical beat signal at the receiver. The information on amplitude /o,thz and phase 
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difference Aip can be separated by phase modulation with a mechanical delay stage or, in 
our case, a fiber stretcher j^. We focus on the amplitude js] 

/0,THz (X P^£,j,, X ^ X (Pl,R.P.,R.)^/^ (2) 

where Pi,pu denote the optical powers of the two lasers at the receiver, and 5'rx and dS'Rx/ dt/ 
are the responsivity of the receiver and its voltage derivative, respectively. Here, we have 
used the linear photocurrent- voltage characteristic of a photomixer (see Fig. 1), thus oc 
-f'THz.Rx X dS'Rx/dt/ with the electric field amplitude -Ethz^x oc PthzRx- Equation ([2]) 
assumes perfect spatial overlap of the two lasers which is exactly the case in our setup where 
the superimposed light is guided in single-mode fibers. However, if the photomixers are 
illuminated by a free beam, spatial overlap might have to be considered. Altogether this 
yields an expression for the amplitude of the terahertz signal. 



-^0,THz oc Tgl/j^pjp X S'Tx(f^bias) ^ (-Pl,Tx-P2,Tx-Pl,Rx-P2,Rx) (3) 



Apparently, /o,thz depends on all of the four optical powers Pjj , which in principle may 
drift independently. A drift of e.g. the optical output power of laser 1 affects Pi_tx and Pi^rx 
but not P2,Tx or P2,rx; whereas an attenuation or mechanical displacement of the beam in 
only the transmitter branch will change Pi,tx and P2,tx- Polarization effects can influence 
all components differently. The polarization plays a role here because of the finger structure 
of our photomixers which more strongly reflects the component of the E field parallel to the 
fingers than the perpendicular one (see e.g. [lo|). 

In order to fully compensate the above mentioned drifts, one would have to measure all 
four optical power components as well as the responsivities (or dS'/df/), which is not easily 
practicable. As a sophisticated alternative, we propose to determine DC photocurrents in 
both photomixers, treating both of them like photodiodes. We will show that this enables 
us to compensate for drifts of the responsivities and of the total optical powers illuminating 
the two photomixers. As a result, the normalized signal is only sensitive to a drift of the 
relative power of the two lasers. 

To the best of our knowledge, biasing a photomixer at the receiver side has not been 







reported lUJ. The transmitter and the unbiased receiver have been discussed above, see 
Eqs. ([1]) and ([2]). Here, we additionally bias the receiver photomixer with a small DC 
voltage t/bias Rjc of typically > 30 mV. Such a small bias does not change the linear behavior 



of the receiver photomixer (see Fig. [T]), thus the shghtly biased photomixer can still be used 
as a detector. For small values of f/bias,Rx we may write the responsivity as S^x{Uuas,Rx) = 
Uhias,RK X dS'pix/dt/. Thus the DC photocurrent Jrx is given by 

-^Rx = f^bias,Rx (-Pl.Rx + P2,Rx) ■ (4) 

Typically, /thz (and /tx) is not a true DC current because the bias voltage at the transmit- 
ter is modulated e.g. in the kHz range in order to apply lock-in detection of Jthz- Therefore, 
the DC photocurrent /pu can be separated from Jthz- 

Even an unbiased photomixer shows a small residual photocurrent due to a photovoltage 
which arises from inhomogeneous illumination Q]. Experimentally, we estimate this offset 
voltage to be of the order of a few millivolt. The offset voltage itself is sensitive to drifts 
of the optical power and thus unsuitable for using the photomixer as a reliable powermeter. 
This task requires a stable voltage, independent of the optical power. Hence, the reliability 
of the normalization is improved if one selects specimen with low offset voltages, employs 
low optical power, and applies a DC bias voltage which is significantly larger than the 
self-induced offset voltage of the photomixer. 

In short, we determine three different currents. The photocurrent /tx cos(ci;mod ^) at the 
transmitter and the terahertz photocurrent Jthz cos(cJmod t) at the receiver both oscillate 
at the frequency Wmod of the bias modulation of the transmitter, whereas Jrx is a true DC 
current due to the DC bias voltage at the receiver. Then, Jtx and Jrx are used to normalize 

-^THz, 

J- _ J -^Tx,0 -^Rx,0 rpl/2 (-Pl,Tx-P2,Tx-Pl,Rx-P2,Rx)"^''^ /^.v 

JTHz.norm - ^THz —j 7 OC J , j,, — ^ TT^ — ^ r [p) 

JTx -(Rx l-r l,Tx + -r 2,TxJl-' l,Rx + -r 2,RxJ 

where /tx,o and /rx,o are constants, e.g. the long-term average values of Jtx and Jrx, respec- 
tively. Note that the normalized photocurrent does not depend on the responsivities 5'tx 
and S'rx anymore. This is valid under the assumption that d^Rx/ df/ is equal for a constant 

1 /2 

f4ias,Rx and for applying a terahertz field -Ethz.Rx cx^ -^thzRx' "which is supported by our 
results described below. 

Let us now examine how this normalization affects the sensitivity to a drift of the optical 
power. The partial derivatives of /thz and /THz.norm with respect to e.g. Pi,tx and P2,Tx are 
given by 

dllYiz _ 1 -^THz <9/tHz _ 1 -^THz 

f^-Pl,Tx 2 Pi Tx dP2,Tx 2 P2 Tx 
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(7) 



Compared to /thz, we find that /THz.norm is less sensitive to a variation of Pi,tx by a factor 
of 

-Pl.Tx — -P2,Tx 
-Pl,Tx + -P2,Tx 

In the hmit of Pi,tx = -f2,Tx; a small drift of Pi,tx is fully compensated. Moreover, 
/xHz.norm ^Iso stays constant for a common attenuation of both lasers, i.e. for power changes 
dPi Tx/-Pi,Tx = dP2,Tx/-P2,Tx- However, full compensation is not achieved if Pi.tx and P2,tx 
are significantly different and, at the same time, change differently. Therefore, equal laser 
powers are preferable for the normalization. An analogous calculation yields the same results 
for the receiver. 

It is instructive to rewrite /xHz.norm in terms of the splitting ratios of the two laser powers 
at the transmitter ttx = Pi,Tx/ (Pi,Tx + P2,Tx) and at the receiver tr^ = Pi,rx/ (Pi,Rx + P2,Rx), 

/THz,norm (X T.^^i, X [rTx(l - T^.f^ [rRx(l - T^^f^ ■ (8) 

The normalized terahertz photocurrent is no longer a function of the total power or of 
the responsivities of the photomixers, but it only depends on the splitting ratios in the 
two branches. In the simple case of equal splitting ratios for transmitter and receiver, 
ttx = '''rx = Eq- (E]) reduces to 

/THz,norm OC T^^^i^ X r(l - r) . (9) 



III. EXPERIMENTAL SETUP 

A sketch of our experimental setup is given in Fig. |2l for details we refer to Refs. js], Q, |9|. 
For both seed lasers, the optical power at a given frequency is actively stabilized to about 
± 1%, but the tapered amplifier and the photomixers are sensitive to temperature drifts. The 
photocurrent at the receiver is preamplified and then digitized to determine the terahertz 
photocurrent Jthz and the DC receiver photocurrent Jrx- The transmitter photocurrent Jtx 
is measured with the help of a 1 kfi resistor in series with the photomixer. 

The signals of both optical frequencies are amplified simultaneously in the tapered ampli- 
fier, and the resulting beam with a particular value of r is sent to both photomixers. Thus 
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the simple case of equal splitting ratios (see Eq. (Q) is applicable to our setup. However, 
small changes of and in the path from the amplifier to the two photomixers cannot 
be excluded. The gain of the tapered amplifier is wavelength dependent due to reflections at 
the surfaces of the amplifier chip. A temperature-induced change of the chip length there- 
fore leads to a drift of the ratio of the laser powers which can typically range from 40:60 
to 50:50, equivalent to r = 0.4 to 0.5. We thus expect that fluctuations of the normalized 
photocurrent /THz,norm are suppressed to below 5 %. 



IV. RESULTS 



As a first example we study the effect of an artificially introduced drift APi of the power 
Pi of one laser (-Pi,tx oc Pi and Pi,rx oc Pi). For this measurement, the tapered amplifier 
was removed from the setup. Note that Pi was not measured directly at the photomixers but 
with a separate photodiode. Experimentally, we started with equal laser powers, r = 0.5, 
and then reduced Pi, keeping the power of the second laser fixed (see Fig. [3]). According to 
Eqs. (|3]) and (|9]), we expect that the terahertz photocurrent Jthz shows the same relative 
drift as Pi, i.e. A/thz/-^thz = APi/Pi, whereas the relative drift of the normalized current 
/THz.norm IS cxpectcd to be much smaller (solid lines in Fig. |3]). This is corroborated by the 
experimental data. However, the measured drift A/THz,norm is slightly larger than expected. 
This is mainly due to the fact that AJthz/-^thz shows additional fluctuations of the order of 
roughly 3%. If we subtract these deviations between A/thz/-^thz and the straight line from 
A/THz,norm/-^THz,norm, then the result is very close to the predicted curve (see crosses in Fig. 
[3]). The remaining difference can be attributed to, e.g., the small power-dependent offset 
voltage at the receiver caused by inhomogeneous illumination. Although the measured data 
do not fully reach the ideal case of the theoretical prediction, the stability of the signal is 
significantly enhanced by the normalization. For instance a reduction of Pi of 20% gives 
rise to only 3% change of the measured value of /THz,norm- 

In the second example, see Fig. HI we changed the splitting ratio r while keeping the total 
laser power Pi + P2 fixed. In this case, /THz.norm follows the predicted behavior r(l — r) 
excellently. 

As a third example, we deployed the complete spectrometer from Fig. [2] and monitored 
the terahertz amplitude at 600 GHz over 2 hours while periodically changing the labora- 
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tory temperature between 22 °C and 24 °C (see Fig. [5]). The measured amplitude varies 
by approximately 15% peak-to-peak, whereas the stability of the normalized amplitude is 
improved to about 4% peak-to-peak (bottom panel). As discussed in Section IIIIl in our 
setup the splitting ratio r drifts between 0.4 and 0.5 as a function of temperature. Accord- 
ing to Eq. OH]), this corresponds to a drift of /THz,norm of 4%, in excellent agreement with 
our data. The DC photocurrents in the transmitter and receiver are correlated to some 
extent (see Fig. [5|, second and third panel), which here is mainly due to drifts of the tapered 
amplifier which is common to both optical paths. But there are also significant differences 
in the fluctuations of Jtx and /rx? substantiating the necessity to measure both currents 
separately. The normalized amplitude is less stable if only Jxx or only is used for the 
normalization. This is evident from the fourth panel of Fig. 0, which shows /THz,norm,Rx 
and /THz,norm,Tx = -^THz (-^Tx,o/-^Tx)^ • In the latter, the factor /rx,o/-^Rx has been replaced by 
Jtx,o//tx in Eq. 

In Fig. |6]we show the 100 % line of the spectrometer in the frequency range from 50 GHz 
to 800 GHz. The measurements were performed directly after switching on the setup. For 
a step size of 1 GHz, a single run took about 10 minutes, and there was a delay of about 10 
minutes between the two runs. The initial change of the temperature gives rise to a drift 
of the 100 % line, which is significantly reduced by the normalization. As discussed above 
for the data of Fig. [5l the normalization is very well suited for temperature-induced drifts. 
However, it does not have any effect on the noise-like features present in the 100 % line. 
These are caused by e.g. the uncertainty in the determination of the amplitude /o,thz from 
the raw data measured as a function of Aip (see Eq. (I3l)) and by small fluctuations of the 
frequency. 

Finally, we varied the receiver bias and measured the noise photocurrent, i.e. the standard 
deviation of the terahertz photocurrent with blocked terahertz beam. We found that the 
noise photocurrent depends on the details of the photomixer device. For some devices, the 
noise is nearly independent of the bias, while other devices show a significant increase of 
the noise photocurrent when a bias of a few 10 mV is applied. Presumably, this difference 
originates in differences in the photomixer resistivity. 

In order to obtain stable and reliable data, it is of course desirable to reduce the fluc- 
tuations of the laboratory temperature in the first place. However, stabilization to much 
better than ±1 K is not a trivial task. Moreover, the stabilization achieved via photocur- 
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rent normalization may be a significant improvement for teraliertz applications requiring 
measurements outside a regular laboratory. As an alternative to the method described here, 
one may consider to monitor the optical power of the two lasers with, e.g., a photodiode. 
However, the discussed normalization via the photocurrents has two main advantages. It 
does not require an extra sensor, and the power is measured directly within the photomixers. 
Therefore, the normalization also compensates drifts of the responsivity of the photomixers, 
e.g. caused by mechanical drifts within the photomixer device. 

V. CONCLUSION 

We have described and demonstrated a normalization scheme for the terahertz photocur- 
rent in a continuous-wave photomixing system. This method is based on measuring the DC 
photocurrents in both the transmitter and receiver photomixers. Consequently, no extra 
sensor is needed and the method can easily be implemented. Any change of the laser power 
illuminating the photomixers can be described as a change of the total laser power in com- 
bination with a changing splitting ratio r of the two laser powers. The normalization fully 
compensates drifts of the total laser power as well as drifts of the responsivities, thus the 
stability of the normalized terahertz signal only depends on the splitting ratio r. We have 
provoked large changes of the terahertz signal by either reducing the power of one laser or by 
changing the laboratory temperature to simulate unstable ambient conditions. In all cases, 
the normalized signal is stable within a few percent. 
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Fig. 1. Current- volt age characteristic of an illuminated photomixer. At low voltages, 
the photocurrent is proportional to the voltage (black points: data, red line: linear fit). 
Consequently, the photomixer converts a terahertz electric field linearly into a photocurrent 
(see Eq. ([3])) and can be used as a detector. Inset: Photocurrent- voltage characteristic for 
larger voltages (see also |12|). 

Fig. 2. Setup of the continuous- wave terahertz spectrometer (from {q]). 

Fig. 3. Drift of the measured terahertz photocurrent Jthz (black) and of the normalized 

photocurrent /THz,norm (red) upon reduction of the power Pi of one of the two lasers. APi = 

corresponds to equal laser powers, i.e., r = 0.5 or Pi = P2. Solid lines: predictions according 

to Eqs. ([3]) and iQ. Full symbols: measured data. Crosses are obtained by subtracting the 

deviation between black symbols and black line from the red symbols. 

Fig. 4. Normalized terahertz photocurrent /THz,nomi vs. splitting ratio of the two laser 

powers, r = Pi /(Pi + P2), for a constant total power Pi + P2 = const. The solid line 

depicts our expectation /THz.norm oc r(l — r), see Eq. /THz,norm is scaled such that 

-^THz,norm('^ = 0.5) = 1. 

Fig. 5. Effect of a periodic variation of the laboratory temperature by ±1 K (top panel) 
on the stability of the measured terahertz signal /thz (black line, bottom panel). The drift 
is strongly suppressed in the normalized photocurrent /THz,nomi (red line, bottom panel). 
Second and third panel: DC photocurrents Jtx and Jpu at transmitter and receiver, respec- 
tively. Fourth panel: Normalized photocurrent if only Jtx (green) or only Jrx (blue) is used 
for the normalization, e.g. /THz,norm,Tx = ^thz (/tx,o/^Tx)^- 

Fig. 6. A 100 % line, i.e. the ratio of the terahertz signals of two consecutive runs, measured 
directly after power up of the system. The data have been smoothed by averaging over 5 
points. 
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FIG. 1: Current-voltage characteristic of an illuminated photomixer. At low voltages, the pho- 
tocurrent is proportional to the voltage (black points: data, red line: linear fit). Consequently, the 
photomixer converts a terahertz electric field linearly into a photocurrent (see Eq. ^) and can 
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used as a detector. Inset: Photocurrent-voltage characteristic for larger voltages (see also 
FiglPhotomixerCurrentVoltage.eps. 
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FIG. 2: Setup of the continuous- wave terahertz spectrometer (from [9(). 

Fig2TIIzSetupFiberStretcher.eps. 
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FIG. 3: Drift of the measured terahertz photocurrent /thz (black) and of the normahzed photocur- 
rent /THz.norm (red) upon reduction of the power Pi of one of the two lasers. APi = corresponds 
to equal laser powers, i.e., r = 0.5 or Pi = P2. Solid lines: predictions according to Eqs. ([3]) and 
0. Full symbols: measured data. Crosses are obtained by subtracting the deviation between 
black symbols and black line from the red symbols. Fig3normPC.eps. 




FIG. 4: Normalized terahertz photocurrent /THz.norm vs. splitting ratio of the two laser powers, 
r = Pi /{Pi + P2), for a constant total power Pi + P2 = const. The solid line depicts our expec- 
tation /THz,norm OC r(l — r), SCC Eq. ([9]). /THz,norm is Scaled such that /THz,norm 

(r = 0.5) = 1. 

Fig4ITHzVsSplittingRatio.eps. 
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FIG. 5: Effect of a periodic variation of the laboratory temperature by itlK (top panel) on the 
stability of the measured terahertz signal /thz (black line, bottom panel). The drift is strongly 
suppressed in the normalized photocurrent /THz,norm (red line, bottom panel). Second and third 
panel: DC photocurrents Itx and /rx at transmitter and receiver, respectively. Fourth panel: 
Normalized photocurrent if only Jtx (green) or only /rx (blue) is used for the normalization, e.g. 
/THz,norm,Tx = -^THz {iTxfl/hxf- Fig5PhotocurrentCorrectionFigPaper.eps. 
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FIG. 6: A 100 % line, i.e. the ratio of the terahertz signals of two consecutive runs, measured 
directly after power up of the system. The data have been smoothed by averaging over 5 points. 
Fig6100percent.eps. 
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